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Modern electronic devices are unthinkable without the well-controlled formation of interfaces at 
heterostructures. These often involve at least one amorphous material. Modeling such interfaces 
poses a significant challenge, since a meaningful result can only be expected by using huge models 
or by drawing from many statistically independent samples. Here we report on the results of high 
throughput calculations for interfaces between crystalline silicon (c-Si) and amorphous silicon nitride 
(a-SisNs.siH), which are omnipresent in commercially available solar cells. The findings reconcile 
only partly understood key features. At the interface, threefold coordinated Si atoms are present. 

These are caused by the structural mismatch between the amorphous and crystalline part. The 
local Fermi level of undoped c-Si lies well below that of a-SiN:H. To align the Fermi levels in the 
device, charge is transferred from the a-SiN:H part to the c-Si part resulting in an abundance of 
positively charged, threefold coordinated Si atoms at the interface. This explains the existence of a 
positive, fixed charge at the interface that repels holes. 


I. MOTIVATION 

Silicon is the most important material for single¬ 
junction solar cells. Although there are other promising 
compounds for photo voltaic devices such as GaAs, GdTe, 
or InP, silicon based cells still maintain about 90 % mar¬ 
ket share m- Silicon has the advantage of being widely 
available, non-toxic, and applicable to many different re¬ 
quirements. The possible applications range from less 
efficient, but cheap, amorphous cells over cost-efficient 
multi-crystalline cells up to highly efficient, but more ex¬ 
pensive, mono-crystalline cells mm- To produce high 
performance cells, optical losses must be minimized and 
recombination centers passivated. Meeting both chal¬ 
lenges, amorphous silicon nitride (a-SiN:H) is commonly 
deposited as anti-reflect ion and passivation layer on top 
of crystalline silicon (c-Si) solar cells with p-type base and 
n-type emitter. Using plasma enhanced chemical vapor 
deposition (PEGVD) [H [5], the deposition parameters 
(i.e. the gas mixture of N 2 /SiH 3 or NH 3 /SiH 3 ) control 
the refractive index as well as the defect concentration 
of a-Si 3 Na,:H so that the reflectivity and the passivation 
can be tuned for optimal cell performance. 

However, even with high quality surface passivation, 
the defect concentration at interfaces is larger than in 
bulk materials, and defect assisted Shockley-Read-Hall 
(SRH) recombination remains a major source of carrier 
losses in solar cells mm- Especially, defect levels in the 
middle of the Si band gap are effective recombination 
centers and carrier traps. They reduce the lifetime of 
electron-hole pairs in c-Si and, therewith, the efficiency 
of the solar cells mm- Earners et al. showed that the 
defect density correlates with what is called the ’’fixed 
charge” [TO]. This fixed charge is positive and located at 
the c-Si/a-SiN:H interface. In contrast to charge neutral 
recombination centers, the fixed charge increases the life¬ 
time of charge carriers in n-type doped c-Si by repelling 
holes from the interface. As undercoordinated Si atoms 
can have negative, neutral, or positive charge (i.e. K“, 


K" , or K+ defects), the fixed charge is typically explained 
by an increased number of K+ defects at the interface 
m- However, it is still unclear how those defects form 
and whether one could also produce c-Si/a-SiN:H inter¬ 
faces with a negative, fixed charge that repels electrons 
and attracts holes. It is therefore obvious that defects are 
important to understand the properties of c-Si/a-SiN:H 
interfaces. Nonetheless, they are very difficult to investi¬ 
gate by experimental methods alone. For instance, life¬ 
time measurements estimate defect concentrations indi¬ 
rectly, and electron spin resonance (ESR) measurements 
can solely detect states occupied by a single electron [12] . 


For this reason, several authors recently examined 
c-Si/a-SiN:H interfaces by using computer simulations 
[T3ll2Q] . For example, Butler et al. performed a topo¬ 
logical analysis and compared the density of structural 
defects for different thicknesses of the transition region 
between c-Si and a-SiN:H [T5l 116]. The study suggests 
an optimal thickness of 2 nm, however, a computationally 
cheap Tersoff potentials was used and, like other studies 
employing force field methods m, their work was limited 
to investigations of structural properties. Specifically, 
their conclusions on electronic properties were based on 
the assumption that all undercoordinated atoms are re¬ 
combination centers. As the electronic properties are im¬ 
portant for most applications, ab initio methods seem to 
be more suitable. Pham et al. generated c-Si/a-SiN:H 
interface structures by combining classical molecular dy¬ 
namics and ab intio methods [20]. Afterwards, they ex¬ 
amined the Go Wo corrected band offsets and the change 
of the dielectric constant perpendicular to the interface. 
Their focus clearly lied on the application of a-SiN:H as 
high K dielectric in MOS-FET devices, since recombina¬ 
tion centers were not examined in their work. 
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II. MODELING SETUP c-Si a-SiN:H 


In the present paper, we explicitly investigate the elec¬ 
tronic properties of gap states at c-Si/a-SisNs.siH inter¬ 
faces with 1 and 11 atni% hydrogen (H). The undoped 
c-Si part is modeled by a hexagonal ^/7 x ^/7 unit cell 
containing 4 double layers of 7 Si atoms stacked in the 
(lll)-direction with a corresponding lattice constant of 
a = 5.43 A. The a-SiN:H part is made up by 59 Si, 69 
N and 1 or 17 H atoms. We have also investigated 4x4 
unit cells with 6 Si double layers of 16 Si atoms finding 
virtually identical results. The respective densities and 
other computational parameters are chosen in accordance 
with our previous study [21]. Since defects are minority 
species which are difficult to describe using force fields, 
our study uses a large scale ab initio molecular dynamics 
approach based on density functional theory (DFT) [2T1 - 
[25| . 90 statistically independent samples are prepared by 
cooling a-SisNs.siH from the melt (2500 — 2700 K) to the 
amorphous state (1500 K) in 37.5 ps. To compensate for 
the lack in system size, we afterwards average over the 
microscopic configurations. Finally, we examine the elec¬ 
tronic gap states and link them to the geometrical prop¬ 
erties of the structures [26]. Before we continue, we note 
that our interfaces are sharp. This has been achieved by 
keeping the c-Si atoms partially fixed by allowing only 
for in plane movements at the interface. Eventually, we 
have removed this restriction at the interface in the final 
relaxation step. 

III. RESULTS 

A. Spatial location of gap states 


In the following, we summarize our most important 
results by investigating the layer resolved number den¬ 
sities and the electronic density of states (DOSs) at the 
c-Si/a-SisNs.siH interface (Fig. [^. The structural mis¬ 
match causes slight distortions in the first double layer of 
c-Si, and more important, long range density fluctuations 
in the a-SiN:H part, with a pronounced increase in the 
N density around 1.7 A (upper panel, red lines) above 
the Si topmost layer. These N atoms always bind to Si 
surface atoms, saturating 60 % of the Si surface dangling 
bonds. The other Si surface atoms bind to Si atoms in 
the amorphous part. This causes the onset of a peak 
in the Si number density. The corresponding peak inte¬ 
grates to 1.3 Si atoms per Si surface atom, and the region 
is unusually Si rich, with very little N content. Analysis 
of the bonding topology suggests that a large number of 
threefold coordinated Si K defects exists in this region, 
similar to, but in much greater number than in bulk a- 
SiN:H [26]. These geometrical defects also cause a large 
number of electronic defect levels located in the band 
gap of c-Si [lower panel, region (in)]. On the c-Si side, 
our analysis at first sight also suggests a large number of 
defects levels [lower panel, region (ii)], however, inspec- 
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FIG. 1: (color online) (a) Number densities and (b) layer 
resolved density of states (DOS) at the c-Si/a-SENs.s:!! in¬ 
terfaces with 1 and 11 atm% H (solid/dashed lines), (a) The 
N number density shows a peak close to the interface inte¬ 
grating to 60 % of the number density in the top most Si 
layer, (b) The considered electronic states lie within ±0.4 eV 
of the Fermi level. The states are found to be predominantly 
localized either at atoms with ideal- or under-coordination 
(red/black arrow). These states are related to the conduction 
band of c-Si or to undercoordinated Si atoms in a-SisNa.siH. 
In some microscopic models, a charge transfer from occu¬ 
pied coordination defects (red isosurface) to the Si conduction 
band (blue isosurface) occurs (green arrow). 


tion of the charge density of these states indicates that 
these are conduction band like c-Si states (red arrow). 
We note that analysis of the large 4x4 unit cell yields 
quantitatively identical results (not shown). 

One central finding is that charge is partially transfered 
from originally occupied K defects in the a-SiN:H to the 
c-Si conduction band. This happens in 7 % of our micro- 
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FIG. 2: (color online) Typical band structure and DOS of a 
c-Si/a-SisNs.siH interface sample in the presence of a single 
defect state. The color coding indicates the degree of lo¬ 
calization and corresponds to the inverse participation ratio 
(red for localized states). Additionally, we show the DOS of 
a bulk undoped c-Si and a defect free bulk a-SisNs.siH sam¬ 
ple. Note that Ep refers to the average Fermi level of the 
interface models. 



FIG. 3: (color online) Two examples for occupied K defects. 
In sample B, the defect level is close to the valence band, 
whereas, in sample G, it is close to the conduction band. In 
sample G, a partial charge transfer into the conduction band 
is observed. Both defect states interact with the valence or 
conduction band states resulting in forbidden crossings. 


scopic models (green arrow). This result, which is sub¬ 
stantiated in detail below, explains the existence of the 
positive, fixed charge at c-Si/a-SiN:H interfaces. When H 
is present, the number density of H is clearly increased in 
the a-SiN:H part at the interface (upper panel, blue line), 
but the N and Si number densities are hardly changed. 
Although H reduces the number of K defects at the in¬ 
terface almost by a factor 2 (lower panel, solid/dashed 
black line), our general observations remain unchanged. 


B. Band structure and defect levels 

We now investigate three exemplary microscopic mod¬ 
els in more detail. The first sample A contains a sin¬ 
gle localized doubly occupied K“ defect in the middle 
of the gap (Fig. [^. In addition to the band structure 
of sample A, we show the DOSs of bulk undoped c-Si 
and a defect free bulk a-SiN:H model structure (left and 
right most panel). These bulk states have been aligned 
at the Is core levels of the Si and the N atoms in the 
middle of the c-Si and a-SiN:H part, respectively. The 
comparison of the bulk DOSs shows that the indirect 
band gap of c-Si (about 0.5 eV) is much smaller than 
the band gap of a-SiN:H (about 2.3 eV), and the conduc¬ 
tion band offset (about 1.0 eV) is larger than the valence 
band offset (about 0.8 eV). That the band gaps are too 
small compared to experiment is a well known DFT ar¬ 
tifact, however, as shown for instance in Ref. im this 
error changes the mid gap levels, which are particularly 
relevant for this study, only very little. The present off¬ 
sets are consistent with photo emission experiments and 
previous simulations on c-Si/c-Si 3 N 4 interfaces pTl \Tf\. 
The band structure furthermore demonstrates that the 
in plane dispersion of the defect level is negligible. This 
indicates that the lateral supercell size is generally suffi¬ 
cient, as there is no interaction between the periodically 
repeated defect images. However, the valence and the 
conduction band of c-Si show a strong dispersion neces¬ 
sitating a dense 8x8 k-point sampling parallel to the 
interface. 

Since defect related states in the amorphous part can 
be situated anywhere within the gap of a-SiN:H [26], they 
often overlap with the c-Si conduction band or valence 
band. To give examples for such situations we show the 
band structure of two samples B and C (Fig.[^. In sam¬ 
ple B, a doubly occupied K“ defect merges into the c-Si 
valence band, and, in sample C, an initially doubly oc¬ 
cupied K“ defect merges into the c-Si conduction band. 
In both samples, the defect states are again strongly lo¬ 
calized and practically fiat. Nevertheless, since small hy¬ 
bridizations between the c-Si and the defect level exist, 
forbidden crossings are observed. The most interesting 
case is sample C. Here the average defect level is slightly 
above the conduction band minimum of c-Si, causing a 
partial transfer of charge from the defect level into the 
c-Si conduction band. This has reduced the occupancy 
(OCC) of the defect from K“ to and relates to the 
main finding mentioned above. The opposite, a partial 
transfer from the c-Si valence band to an originally empty 
K+ defect, is never observed. 


C. Electronic states in the band gap 

To obtain a statistically meaningful result, we now ex¬ 
plore the electronic density of states evaluated close to 
the interface and averaged over all 90 samples (Fig. [^. 
On the c-Si side, we observe a slight band tailing intro- 
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FIG. 4: (color online) Local electronic density of states eval¬ 
uated for the c-Si part and in the a-SiN:H part (red/black 
lines) of the c-Si/a-SisNs.siH interfaces with 1 and 11 atm% 
H (solid/dashed lines). The mid gap level of bulk c-Si and 
the Fermi level of a-SisNa.sHo.s are shown by vertical lines in 
the corresponding colors. 


duced by local disorder and distorted Si bonds (red lines). 
However, the majority of the gap states are caused by 
undercoordinated Si atoms at the a-SiN:H side of the in¬ 
terface (black lines). The most important observation is 
that the minimum of the DOS of c-Si is shifted towards 
lower energy values (left) in comparison to the minimum 
of the DOS of a-SiN:H. Likewise, the mid gap energy level 
of bulk c-Si (red vertical line)— evaluated by taking the 
mean value of the valence band maximum and conduction 
band minimum of bulk Si — lies about 0.16 eV below the 
average Fermi level of bulk a-SiN:H (black vertical line). 
Impaired by the issue of too small DFT band gaps, we 
believe that we observe only the lower bound of this off¬ 
set here. Note that we have again taken the core levels of 
Si and N as reference for the alignment, and the general 
observations apply independently of the H concentration. 
As for our previous studies on bulk a-SiN:H [211 ES], H 
mainly reduces the number of defects, although here it 
also slightly reduces the misalignment of the Fermi levels 
from 0.24 eV to 0.16 eV. 

The offset of the mid gap level of Si and the Fermi level 
of a-SiN:H has important consequences. In a macroscopic 
sample, the Fermi levels must align and, in a real de¬ 
vice, the alignment would occur via charge transfer from 
donors, here doubly occupied K“ defects in a-SiN:H, to 
acceptor levels in c-Si as illustrated in Fig. [^ In our 
microscopic realizations this is not strictly so, since the 
Fermi level of the c-Si part can be located anywhere be¬ 
tween the conduction and valence band: in the absence 
of acceptors and donors, the c-Si part does not pin the 
Fermi level mid gap. Only if an occupied K“ defect is lo- 
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FIG. 5: (color online) Illustration of the predicted band 

bending at the c-Si/a-SiN:H interface. Before bringing both 
subsystems in contact, the mid gap level of c-Si lies below 
the Fermi level of a-SiN:H (offset between red and blue line 
at the interface). In a-SiN:H, nearly all defect levels below 
and above the Fermi level are doubly occupied K“ defects 
and unoccupied defects, respectively. Only a few singly 
occupied K° defects are located at the Fermi level. To align 
the Fermi levels, K“ and K° defects donate charge to shallow 
acceptor levels in c-Si (red area). They are thereby converted 
to K° and defects (blue area). 

cated above the conduction band of c-Si, a partial charge 
transfer to the c-Si part is observed, like in sample C. This 
partial charge transfer shifts the c-Si states upwards and 
the a-SiN:H states downwards. This line of thought in¬ 
dicates that it is not quite straightforward to determine 
highly accurate conduction band and valence band off¬ 
sets between amorphous model structures and crystalline 
samples by using small model structures, even though 
we average over many microscopic models. That the mid 
gap level of c-Si lies well below the Fermi level of a-SiN:H 
is, however, a robust feature corroborated by three inde¬ 
pendent observations: (i) the observation of microscopic 
models with charge transfer to c-Si, (ii) the visual mis¬ 
alignment of the minima in the electronic DOSs, and (iii) 
the misalignment of the calculated c-Si mid gap level and 
the Fermi level of a-SiN:H. 


IV. CONCLUSIONS 

In summary, we conclude that the number of underco¬ 
ordinated Si atoms (K defects) is significantly increased 
at c-Si/a-SiN:H interfaces compared to bulk a-SiN:H 
(compare Fig. Bb)). This is a result of the structural 
mismatch between c-Si and a-SiN:H. The number of co- 
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ordination defects and related electronic defect states is 
reduced by about a factor 2 in our simulations when H 
is present. Hence, hydrogen clearly cures geometrical de¬ 
fects at the interface, similar to bulk a-SiN:H. We predict 
a positive, fixed charge at c-Si/a-SiN interfaces related 
to a surplus of K+ defects on the a-SiN:H side screened 
by negatively charged acceptor states on the c-Si side 
(Fig.§. In principle, this surplus is related to the mid 
gap level of bulk c-Si lying at least 0.16 eV below the 
Fermi level of bulk a-SiN:H. Consequently, electrons must 
be transfered from originally doubly occupied K“ defects 
on the a-SiN:H side, to acceptor levels on the c-Si side, to 
align the Fermi levels in the heterostructure. Solving the 
Poisson equation for this standard device problem results 
in a band bending at the interface (indicated by black 
lines in Fig.|^. Clearly such a band bending will attract 
electrons to the interface aiding the electron removal by 
metallic leads, which are etched into the a-SiN:H. Holes, 
on the other hand, are repelled making an effective sep¬ 
aration of both carriers possible. Our present study also 
indicates that the number of K defects and the amount 
of positive charge at the interface are interrelated quan¬ 
tities. If the defect density is very low, there will be sim¬ 
ply no defects that can donate charge to the c-Si part. 


This will reduce the positive interface charge and con¬ 
sequently the carrier separation. If the defect density 
is high, the fixed charge density will be high, but more 
recombination centers will be present as well. These re¬ 
sults are exactly in line with recent experimental studies 
[IQIIII]. Based on our present results, we suppose that it 
will be hard to achieve a negative, fixed charge by using 
the standard a-SiN:H passivation. The mid gap level of 
c-Si is well below the Fermi level of a-SiN:H, and only 
strong doping might change this. Moreover, it is impos¬ 
sible to reduce the density of carrier traps (K“, K^, and 
K+ defects) without changing the field effect passivation 
of fixed charges, as the species responsible for the fixed 
charges and K defects are identical. 
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